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1. Introduction. 

This paper studies the propagation of electromagnetic waves through perturbed periodic 
media with period /i ^ lu We treat the resonant case where the length scale of the 
periodic structure is comparable to the wavelength. The observation time t satisfies t ~ 



""^A word on units. The Maxwell equations in vacuum have permittivities e, /i and speed of light 
c = l/y/ejl. Since the speed of light is ^ 1 in KMS or CGS units, ejjL is small in those units. We 
perform an asymptotic analysis as /i — >■ 0. Denote by At the unit of time. No matter what the units 
one has h ^ At/{e, ^} in this limit so there is scale separation no matter what are the values of e and 
jjL. Nevertheless it is wise to, and we choose to, work in units with cAi comparable to 1. For example 
centimeters for length and At ^ 1 equal to the number of seconds that it takes light in vacuum to 
traverse one centimeter. In those units cAi = 1 and one expects that the constants in our error bounds 
will not be very large. 
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\/h. This is the diffractive time scale where standard Maxwell equations (without periodic 
structure) are approximated by Schrodinger's equation (see [12] )• Wavelengths that are 
short compared to the period are short compared to the scale on which the coefficients 
vary. This is the domain of validity of standard geometric optics (see [12], [7] for the 
diffractive case). Wavelengths long compared to the period are analysed by standard 
homogenisation [8]. The interest, both mathematical and scientific, of the resonant scaling 
is that the speeds of propagation and diffractive effects for wave packets are given by the 
Bloch dispersion relation of the periodic medium and not by the symbol of the hyperbolic 
operator or its hyperbolic homogenisation. The propagation speeds can be radically 
different from those of the original equations. The new speeds must not violate the 
finite speed of the original equations (see Theorem 17.21 for a proof of this upper bound) 
but can be much smaller. This is the basis for strategies to slow light ([S], [TB], [B], [22] )• 
That in turn is one of the proposed design elements of the all optical computer. Another 
domain of application is photonic crystal fibers constructed with periodicity in crossection 
(see [27], [H], [21]). The last three examples are modeled by Maxwell's equations that 
are the subject of the current article. 

In our papers [1], [2| we study scalar wave equations. In most cases. Maxwell's equations 
cannot be reduced to scalar equations. The most interesting applications require the 
methods of the present paper. For constant scalar permittivities the Maxwell equations 
can be reduced to the scalar wave equation. For constant e, with e a three by three matrix 
with three distinct positive eigenvalues, it has been known since the time of Hamilton (see 
|15] . |11] page 610), that the characteristic polynomial is of the form t'^Q{t,C,) with Q 
an irreducible quartic polynomial. The characteristic variety is conic with nontrivial 
singular points. Q does not factor as the product of two quadratics in which case the 
variety would be the union of two (double) cones with elliptic cross-section. For variable, 
even scalar, permittivities, the standard derivation of second order equations by taking 
time derivatives works but leads to a system of second order equations for E coupled 
through lower order terms. 

Maxwell's dynamic equations for unknown E{t,x), B{t, x) G X read 



whose infinitesimal generator is not elliptic. Indeed if e, /i depends only on x, then 
(Va;0, Vx4') is a stationary solution for any 0(x),'?/'(x) G C^(M^). For any ball there 
is an infinite dimensional set of such solutions supported in the ball. This is one of the 
main differences of the current article with the earlier ones on scalar wave equations. The 
failure of ellipticity is compensated by the fact that the set of physically relevant solutions 
have additional strong control on their divergence. Those solutions satisfy semiclassical 
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ellipticity estimates (see Theorem 13. 2p . A second principal difference with the earher 
work is that for systems it is not uncommon to have Bloch eigenvalues of multiplicity 
greater than one. In that case the Schrodinger equations of diffractive geometric optics 
are systems and we treat that possibility. A third difference with the scalar case is a 
simplification. For the scalar case, estimates for gradients of the error in the approximate 
solution were straight forward but there was a strikingly difficult argument to estimate 
the undifferentiated error. In the case of systems, the natural energy estimate is an 
estimate. Derivative estimates are proved from Li^ estimates for the time derivatives and 
the divergence. The remaining derivatives are estimated by an ellipticity argument. The 
difficult argument from the scalar case is not required. 

There are two problems closely related to the ones we analyse. The first is the treatment 
of first order elliptic systems. The second is the treatment of the Maxwell system on the 
time scales of geometric optics. The second problem is solved en passant in ^ We have 
chosen to skip the first and jump directly to the Maxwell's equations. The methods that 
suffice for Maxwell yield the elliptic case directly. 

The semiclassical estimates of the present article permit a strengthening of the earlier 
paper [2], where derivatives of order < 1 were estimated. The new method yields estimates 
for derivatives of all orders. 

Our three articles use corrector terms in asymptotic expansions. Article |I] correc- 
tors were constructed by an ad hoc method and used in test functions to study weak 
convergence. In [2] we introduced a a general strategy yielding accurate expansions. 

The /i-dependent Maxwell equations are written in the form P'^lE^, B^) = with 



For the diffractive scaling the perturbations satisfy the following hypothesis where eo{x/h) 
and fiQ^x/h) are the permittivities of the unperturbed periodic structure at scale h. The 
6x6 matrix valued function serves for example to model dissipative effects such as 
Ohm's law (see §H]). 

Notation. For two vectors (e, 6) G C'^ x C^, their Hermitian inner product is denoted by 
(e, b) while (e, b) denotes the ordered pair in C'^ x C^. The cross product in is denoted by 
A. The set ofhnear maps (homeomorphisms) on a vector space K is denoted by B.om{K). 
For any a = (oq, ai, 02, c^s) ^ N'^ the notation 9"^0(t, x) means (IlLi ^x'i) 4>{t, x). 

The Maxwell equations (II. 2p are a symmetric first-order hyperbolic system for = 



(1.2) 
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with symmetric matrix coefficients Aj for j > 1 defined in fl8.3p . f l8.4p and 

Hypothesis 1.1. Diffractive time scale hypothesis. The coefficients in (11.21) are 

given by 

e'^{t,x) = eQ{x/h) + h'^ei{t,x,x/h) , 
H^{t,x) = fio{x/h) + h'^fii{t,x,x/h) , 
M^{t,x) = hM{t,x,x/h) . 

The matrix valued functions eo{y) , fio{y) ^ C°°(T^) are symmetric and positive defi- 
nite and for all a, 5"^ ^{ei, /ii, M}(t, x, G L°°(]R-'^+^ x T^) with denoting the three- 
dimensional torus (R/27rZ)^. Define 

Definition 1.2. The space of Lf^^{M.^) periodic functions of period 2tt is denoted L"^ (T^) . 
The functions {E,B) of L'^{T^) with values in M.^ x M'^ are normed by 

[ (\E\^ + \B\^)dx. 

When normed by the equivalent expression natural in the context of Maxwell's equations 

[ ({E,eoE) + (5,/xo5)) dx 

it is denoted ^^^^^^(T^). 

Definition 1.3. For 6 G [0, Ip a function g{x) on is ^-periodic when x — > e~^^'^ g{x) 
is periodic in x with period 27i. The parameter 9 is called the Bloch frequency. The set 
of Lf^^{E?) 9-periodic functions is denoted L^(Tg). With alternate norm as in Definition 
ITEit IS denoted Ll^^^^{Tl). 

Our solutions are amplitude modulated Bloch plane waves. The plane waves are 
^-periodic solutions of the unperturbed periodic Maxwell equations of the form u = 
e^*{E[x), B{x)). Equivalently E,B are solutions of the spectral problem 

We recall in §21 that the spectrum at fixed 6 consists of {0} with infinite multiplicity 
and a discrete set of purely imaginary eigenvalues A = iu{6) of finite multiplicity. We 
label the nonzero eigenvalues according to their distance from the origin and repeat them 
according to their multiplicity 

(1.5) ••• < 10-2(0) < uj-i{9) < uJo{9) = < uji{9) < 102(9) < 



We work near an eigenvalue of constant multiplicity. 

Hypothesis 1.4. (^Constant multiplicity hypothesis.^ Fix6_ 7^ 0, n G Z*, and denote 
by K the multiplicity, 

Assume that there are 5i > 0, ^2 > and real analytic u}{6) defined on {\6 — 6] < 61} so 
that for each \0 — 6] < 61 the only point of the spectrum in {X : \X — iUn{0.)\ < ^2} is iuj{6) 
with multiplicity k. 

The hypothesis is automatically satisfied when k = 1. An example with k > 1 is e and fi 
constant and scalar where every eigenvalue is of constant multiplicity two (see Example 
below) . 



Definition 1.5. When the constant multiplicity hypothesis \1.4\ is satisfied the group 
velocity is defined by 

V := -Veuj{e). 



Define 

fxoiy)) \-it9 + dy)A 



(1.6) h{co,9,y,dy) := iu 



L with domain equal to the periodic functions in H°°(Ty) := ns>o-ff^(T ) is formally 
antiself adjoint on L^(T'^; dy). The method of proof of Proposition 16.31 shows that the 
closure has domain equal to the v G L^(T^) so that hv G L^(T^) and that the closure is 
antiselfadjoint. 

Definition 1.6. Denote by U the projection operator onto IK := keT'L{u{6),6_,y,dy) along 
the image o/L. 11 is orthogonal with respect to the scalar product o/L^(T^ ; dy) and not 
with respect to the scalar product of L^^^ ^^^(T'^). 

Our wave packets have group velocity V and travel for long times. They see the coeffi- 
cients on group lines for long times. The averages of the coefficients along such long rays 
are important. A particular combination enters in the asymptotic description. 

Definition 1.7. For each t,x define the linear map '~f{t,x) G Hom(]K) by 

(1.7) 7(t,x) := (nA°(t,x)n)"^ n(iu;Aj(t,x) + M(t,x))n. 

We assume that the ray averages 

1 r 

(1.8) lim — / 'y{t,x + Vt)dt := 7(0;), exist uniformly in x G M^. 

We make a fairly weak assumption asserting that this limit is attained at an algebraic 
rate. 
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Definition 1.8. The function satisfies the ray average hypothesis when (11. 8p holds 
and there is a < /3 < 1 so that for all a E N x the solution ga(t, x) of 

(^dt + VA)9a = dr^,{^{t,x)-^{x-Vt)), g^{0,x) = 

satisfies {t)-^g^ e L°^([0, oo[xM3) where (t) := (1 + ^2)1/2, 

This hypothesis, introduced in [2], is discussed in §9.1[ Our main theorem gives an 
approximate solution and an error estimate. In the theorem, T is a new variable, a slow 
time. In the approximate solution it is replaced by ht. In addition there is a K valued 
function, wqITjX). Abusing notation, the value at (T, x) is a function of y denoted 
wo{T,x,y). 

Theorem 1.9. Assume that 7 satisfies the ray average hypothesis with parameter < 
P < 1. For f e S{R^;K) define Wq := Wo{T,x - Vt), where Wq e C"^(M; ^(M^;^)) zs 
the unique solution of the initial value problem for Schrddinger's equation 

(1.9) (^dr + ^td^Mdx,d,)+^ix)')wo = 0, ^o(0, x) = /(x). 

Define a family of approximate solutions 

v''{t,x) := e'("(^)*+^-")/'^ w;o(/it,t,x,x//i) 
and let v!^ denote the exact solution of P^u^ = with u^\t=o = ll^\t=o- Then 
VT > 0, a G N^ 3C(a,T), V/i G]0,1[, sup \\{x , hdt,^)''{u''-v!')\\^, , < C h^~^ , 

with < (3 < 1 from the ray average hypothesis of Definition \l.Si 

Remark 1.10. i. The operator UAqU is positive definite on K. When M = 0, mulitplying 
the Schrddinger equation (1.9) by IIAqII shows that J {HAqUwq , wq) dx is conserved, ii. 
The principal part of equation (11.91) is scalar. Coupling occurs through 7 (see Remark 
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2. L^(R^) ESTIMATES 

Suppose given real symmetric permittivities 

e'^(t,x) , G C°°(M^+^x]0,l[;,; Hom(M3)) 

satisfying the positivity constraint 

(2.1) Vt,x,/i < cJ < e'^(t,x) , /(t,x) < CI. 

Consider the dynamic Maxwell equations in a medium that varies on scale < h << 1 

(2.2) dt{e\t,x)E^) = curl5^ dt{fj.\t, x)B^) = -curlE\ 
This is a symmetric hyperbolic system. The energy identity for solutions is 

(2.3) dt f {E^,e^E^) + {B^,^i^B^) dx = - [ {E^,dte^E^) + {B^ , dt^i^ B^) dx . 

Large time derivatives of e^,fi^ can lead to rapid growth of energy. On the other hand 
if dte^^dt^^ are bounded (resp. 0{h)) one has uniform estimates for the norm for 
t = 0{1) (resp. t = 0{l/h)). Corresponding estimates for derivatives is subtle because 
the coefficients are rapidly varying in space. Our strategy is to derive estimates for 
time derivatives and for divE''^, divi?'^. Then estimate spatial derivatives using an elliptic 
estimate from the next section. 

3. COERCIVITY 

The generator of the dynamic equations (11. 2p is not elliptic. However, the special 
structure of Maxwell's equations yields supplementary bounds on divE and divB. The 
over determined system consisting of the generator together with divergence is elliptic. 
For problems with coefficients oscillating in x on scale h, estimates in semiclassical Sobolev 
spaces are natural. 



Definition 3.1. The semiclassical Sobolev norm if™(M ) denoted \\ ■ ||_f/™(R3) defined by 

1/2 



J2 j \{hd^Tu\^ dx 



\a\<m " 

For a function u{t,x), integer m > and h g]0,1[ define the semiclassical norm with 
derivatives in space and time, 

(3-1) imt, ■■= E ii(^^*>^)"«wiil^(M3)- 

\a\<m 

Denote by C^(Mi+3) the set of families {w'' : < h < 1} C C'=(Mi+3) so that 
(3-2) sup II(^^*.^)°^^Hl-(mi+3) < 

o<h<i ,~r, ^ ' 

\a\<k 
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The same notation is used for functions valued in any finite dimensional real or complex 
vector space. The left hand side of (13.21) serves as norm making Cj^ a Banach space. The 
Frechet space is defined as H/cC^. 

Theorem 3.2. For e'^{t,x), fi'^{t,x) G C^(M^^^) and m G N there is a constant C{m) so 
that for all E,B e H^'iR^), t G M, and h > 

II^WL-(M3) < C{m) (^\\hcuilE{t)\\^^,^^^,^^hdw{e''{t,x)E{t))\\^rn-i^^^^^ , 

II^WL-(R3) < C{m) (^\\hcni\B{t)\\j^rn-i^^3)M\hdiv{fx^{t,x)B{t))\\^rn-i^ . 

Proof. The inequahties for E and B are identical so it suffices to consider E. The variable 
t is simply a parameter so it suffices to consider t fixed. The key estimate is the following. 

Lemma 3.3. If t{x) is a symmetric matrix valued function so that e>cl>0 for all x, 
and d"e G for all a G N'^, then for each m G N there is a constant C{m) so that 

for all E G H'^iM.^) 

(3.3) ||^||^„(K3) < C(m) (||curl^||^_,(^3^ + ||div(e(x)E)||^_,(j^3^ + ||E||^_,(^3^) . 
Proof. The integrand in 

/ |div(e(x)E)|2 + |curl^|2 dx 

is a quadratic form in E and its first derivatives. The terms quadratic in the derivatives 
of E have the form 

^ {aij{x)diE , djE) 

l<jj<3 

with uniquely determined real matrix valued functions aij{x) with aj^i equal to the trans- 
pose of Qij. Introduce the the symbol The definition implies that for any 
X, e and ^ in 

(3.4) 5^(a,,(x)e.0e,e> = |(e,e(x)e)|2 ^ ^^^^^2^ 

Choose < c < 1 so that for all x, e{x) > cl. The lemma is a consequence of Garding's 
inequality once we prove that for all real x, e, ^ 

(3.5) l(e,e(x)e)p + \^ A e\' > c\^\'\e\y2. 

By homogeneity it suffices to consider |^| = 1. 

Decompose e = ey + e^ into parts parallel and perpendicular to ^. The definition of c 
yields 

(3.6) |e^p<|e|V2 =^ | , e(a;)e||) p > c|e,||2 = c{\e\^ - \e^\^) > c|e|V2. 
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On the other hand since c < 1 

(3.7) |exP>|e|V2 =^ 1^ A ej^ = \e^\^ > \e\^/2 > c\e\^/2. 

Estimates (13.61) and (13. 7p imply (13. 5p completing the proof of the lemma. □ 
Scaling shows that (13. 3p implies the h dependent coercivity of Theorem 13.21 □ 

4. Stability 

Theorem 4.1. Let P'^ , he the operator ([OD with , dte^ , ^i^ , dt^i^ , G C^(R^+3). // 
f,g & H°°(M.^) := DgH'^eR^^) then there is a unique family of solutions = (i?^,i?^) G 
C°°(R; if~(R3)) to P^u^ = 0, with u^{0) = {f,g). For each m E N there is a constant 
c(m) so that with 



L°°(Ri+3) 



C{m,h) := \\{hdt,.ndte\dt^\M'')\ 

\a\<m 

one has 

(4.1) \\u^it)\\ , < c(m)e'=('")^('"''')* ||m^(0)|| 

^ ' II ^ ' \\m,h — ^ ' II ^ ' \\m,h 

Remark 4.2. If = 0, one has dt{div{e^it,x) E'')) = and dt{diY{fi^{t,x) B'')) = 0. 
In particular one has 

(4.2) div{e^{t,x)E'') = 0, div(/(t, x) 5^) = 0, 
as soon as these identities hold at t = 0. 

Remark 4.3. For the analysis on the difjractive scale, the theorem is applied with C{m, h) - 
0{h) as h 0. In that case one has uniform bounds fort = 0{l/h). 



Remark 4.4. This careful accounting of derivatives in Theorem 4jJ_ heart of 

extenting the results of this paper to equations whose coefficients are only finitely differ- 
entiable. 

Proof. The existence for fixed h is classical. The estimate for m = follows from Gron- 
wall's inequality together with 

dt I ((e^ , E^) + (/ , B^)) dx = - f {{dte^ E^ , E^) + {dt^i^ B^ , B^)) dx 

+ 2 / {u^,M^u^)dx 

<C(||9ie"(t),9i/(t)|Uoc(R3) + ||M^(t)|Uoc(R3))||(E\5'^)(t)|U2(K3). 
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The proof for arbitrary m is by induction. Assume the estimate proved for m. Denote 
by K'^(t,s) the evolution operator associated to the equation P'^w = 0. Precisely, for a 
distribution g G V'{R^), w{t) = K^{t, s)g is the unique solution of the Cauchy problem 

P^w = 0, wl^^ = g. 

The inductive hypothesis is a bound 

\\K^it,s)g\U^,, < c(m)e^('")^('"'^)l*-^l||(7||H,"^(M3). 
The Duhamel relation 

v^^it) = ir^(t,0)t;^(0) + [ K''{t,s) P^\v%s) ds 

Jo 

then yields the estimate 



The key point is that one cannot simply apply the estimate at level m to := hdu^. 
Doing so yields 



Write P^{hd) = {hd)P^ + [P'^,hd]. The first term vanishes when applied to u. The 
commutator [P^, hd] is 

If d were a derivative with respect to x then dxie'^yfi^} ~ l/h leading to an unac- 
ceptably large contribution. Instead of estimating all derivatives we estimate only dt, 
diy {e^ E^, fi^B^}, and curl. The time derivatives of e^,fi^ are bounded. The divergence 
and curl play a special role in Maxwell's equations. The remaining spatial derivatives are 
recovered using coercivity. @ 
For dt compute 

[P^hdt] = + 

Estimate 

\\P\hdtu'')is)\\m,h < Cim+l,h)\\u\s)\\m+l,h. 

The Duhamel estimate yields 
(4.3) 

\\hdtu\t)\U,f, < c(m)fe^('")^(^'"''^)* 11^/^(0) |U+i,,+ / e^^^'")^^"^'^)^*-^) C{m+1, h)\\u\s)\\m+i,h ds 



^The use of dt, div, curl is reminiscent of the use of dt + 'vdx, div, curl , and tangential derivatives for 
the inviscid compressible Euler equations in 24 . 
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Write the Maxwell equations as 

{hcmlB^ , -hcmlE'') = hdt{e^ , B^) + h . 

Therefore 

\\hcni\E^{t),hcni\B^{t)\[^^ < (^r.h.s. of (gS]) + \\hu^{t)\\m,h)c{m, h) . 
Use the fundamental theorem of calculus to estimate 

\\hu''{t) \\m,h< \\hu''{0) \\m,h + [ \\hdtu''{s)\\m,h ds . 



Wasting a derivative in the first term yields 

< ||m^(0) \[m+l,h+ [ \\u^\s)\\m+l,h ds. 

Jo 

Combining the last four assertions yields 

(4.4) curl , /i curl 11^^ < r.h.s. of (gSD • 

Next estimate \\hdive^E'^, hdiv fj.'^B'^\\m,h- Write in block form with 3x3 blocks 

^ - [m,\ m^J ■ 

Therefore 

dt{diY{e''E'')) = diY{e^ E^)t = -diviM^^E^ + M^^B'') . 
Integrate to find 

\\hdw{e''E'')it)\\m,h < ||/idiv(e^£;'^)(0)|U,, + f \\M^,E\s) + M^,B\s)\\m,h ds 

Jo 

< C(m,/i)||n^O)|U+i,, + / C{m,h)\\u\s)\\^,hds. 



Performing the analogous estimate for hdiv{^^B'^) and replacing || \\m,h on the right by 
the larger || ||m+i,/i yields 

(4.5) \\h diY{e^ E^){t) , hdiY{fi''B^)\\m,h < r.h.s. of (gSD • 

Combining (14. 3p . (14.41) and (14. 5|) . the inductive hypothesis, and the coercivity estimate 
from Theorem 13.21 yields 
(4.6) 

lk'(^)|^n+l,/^<c(m)f||M'*(0)lU+l,,+ / e^('")^("^''^)(*-^)C(m + l,/i)||M'^(s)l|„+i,,(is). 



The theorem follows from Gronwall's inequality. □ 
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5. Stationary solutions 
When e and fi depend only on x, there is a conserved norm for the solution of (11.11) . 

(5.1) dt / {E,e{x)E) + {B , fi{x) B) dx = 0. 

Equations (II. ip have an infinite dimensional space of stationary solutions. The set of func- 
tions satisfying dive(a;)-E = div iJ,{x)B = is invariant and orthogonal in the conserved 
scalar product to the stationary solutions. 

To prove this assertion we use the Fourier Transform. For any u G -/^^(R'^), we have 
u(^) = (27r)-3/2 f e-'""-^ uix) dx, m(x) = (27r)-3/2 f e'^-^u{i)di. 
The Fourier Transform is unitary on L^(M'^, dx). 

Theorem 5.1. When e{x),fi{x) depend only on x, u{x) = {E{x), B{x)) G L^(M^) is a 
stationary solution of the dynamic Maxwell equations (11.11) if and only if curl E = curlS = 
0. The orthogonal complement of these data in L^^(M^) normed by J {eE, E) + (/i-B, B) dx 
is invariant under the flow and consists exactly of the solutions satisfying (14. 2 p . The fields 
(grad(/>, grad^/)) with (f),ip & H^{M.^) are L^(R^) -dense in the stationary solutions. 

Proof. The first assertion is obvious. 

The orthogonal complement of the stationary solutions is invariant because the evolu- 
tion is unitary. 

Next prove the density of gradients. The field E satisfies curlE' = if and only if 
^ A -E(0 = O5 that is the Fourier Transform E{^) is parallel to ^ for almost all ^. For 
^ 7^ E = ^f{^) uniquely defines the scalar valued /(O- Choose a smooth cutoff function 
< X ^ 1 vanishing on a neighborhood of = and identically equal to 1 outside 
a compact set. Then E is the L^-limit as n tends to infinity of the field with Fourier 
transform equal to xi^O C/(0- Define 0" = x(^0/ so G iJ^(]R^) and grad0" E in 
L\ 

Using the density of the Schwartz space 5(M^) in if^(M^), shows that a vector E is 
orthogonal to the stationary states if and only if 

/ (eE,grad0) dx = 0, for all G S{R^) . 

This is the definition of div{eE) = in the sense of tempered distributions. □ 

6. The theory of Floquet and Bloch 

6.1. Bloch Transform. We recall the essentials of the method of Floquet and Bloch (see 
for example [ISj, [91, HO], [26], [30], [8]). 
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de 



Write each ^ G M'^ as n + ^^ with uniquely determined n G Z'^ and 6 G [0, Ip. Expressing 
u{x) in terms of its Fourier transform, u{C), yields 

(6.1) u{x) = {27t)-^/^ [ e'^-" f V e'"-"u(^ + n)) 

The parentheses enclose the Fourier series expansion of a function periodic in x with period 
27r. Considered as a function of x, the integrand is ^-periodic in the sense of Definition 
II. 3[ Identity (16.11) decomposes L^(M^) as the direct integral over 9 of the Hilbert spaces 
of ^-periodic functions belonging to L^q^(R^). 

Proposition 6.1. The map that associates to u its Bloch wave expansion uq 

(6.2) L2(R3)3^ ^ u^^x) := {2ti)~^'^ e'^-"" e'''-''u{e + ti)^ G /.^(T^) 

yields a unitary decomposition of L'^{M.^) as the direct integral over 6 G [0,lp c»/L^(Tg). 
The inverse is given by 

|2 / ||„. ||2 



u{x) = / U0{x) dO with ||m||^2(|53) = / ||'"e||L2(x3) dO . 
J[o,iP ■>'[o,ip 

The map u y-^ e~^^-^ue is a unitary map L^(M^) — )■ -L^([0, Ip ; L^(T^)). 

Remark 6.2. The partial derivatives and 2t^ -translates of 6-periodic functions are 9- 
periodic and the product of a 6-periodic function by a 2'K-periodic function is 6-periodic. 
Therefore partial differential operators with 27i-periodic coefficients map 6-periodic func- 
tions to themselves. The Bloch decomposition reduces these operators. 

6.2. Maxwell's equations. The method of Floquet-Bloch applies to Maxwell's equa- 
tions (see for example ^ and [28]). The delicate part for us is the infinite dimensional 
kernel and the degeneration of the coercivity estimates as /i — )■ 0. From here to the end 
of this section the Bloch strategy is used to analyse Maxwell's equations 

in the case of periodic e{x) and fi{x). It suffices to analyse its action as a map on L^(T^) 

Proposition 6.3. If Aj are symmetric matrices then the operator L = Ajdj satisfies 

(6.4) \/u,v e L\Tl)nC°^ {Lu,v) = -{u,Lv). 

Denote by L the closure of the operator so defined and by L* the Hilbert space adjoint. 
Then L is antis elf adjoint with domain equal to the set of u E i^^(Tg) so that J^j ^j^j'^ ^ 
L^(Tg) in the sense of distributions. 
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Proof. To prove ( 16 ■4p . write u = e^^'^u and similarly v with periodic u,v. Then, 
{Lu,v) = [ {AjdjU,v)dx= I {Ajdj{e'^-''u) , e'^-''v) dx 

{e'^-''Aj{dj + iej)u , e''^-^v) dx = I {Aj{dj + iej)u , v) dx 

[0,27rP i[0,27rP 



{u , Aj{dj + iOj)v) dx , 

[0,27rp 

the last step by integration by parts and periodic boundary conditions. 

Identity fl6.4p implies that L* D — L in the sense that the left hand side is an extension 
of the right. The definition of distribution derivative implies that L*u = / G L^(Tg) if 
and only if —{^j Ajdj)u = / in the sense of distributions. 

In that case denote by Js a standard mollifier. Since 6'-periodic functions are invariant 
under translations, Jsu G C°° fl L^(Tg). Since Js commutes with Ajdj so := Jsu 
satisfies —Lu^ = f. Passing to the limit shows that u belongs to the domain of L and 
Lu = -/. Thus L* C -L. □ 

The spaces L^(Tg) depend on 6. The following proposition allows one to apply standard 
results in perturbation theory. 

Proposition 6.4. The unitary map L'^(T^) 3 v ^-^ g^e.x^ ^ L^iJ'e) intertwines the 
operators 

d^A\ 1 (d^ + i0)A 

' and ' 



-d^A J ' \-{d^ + te)A 

The former acts on L'^iT'g) and the latter on the 9 independent space L^{T^). The latter 
family of operators depends analytically on 6. 

Proof. The unitary map commutes with multiplication operators and intertwines the an- 
tiselfadjoint operators dj on L'^iTl) and dj + iOj on L^(T'^). This yields the desired 
result. □ 

The straight forward proof of the next result is omitted. 

Proposition 6.5. A function u{x) = {E{x), B{x)) G L^(R^) is a stationary solution of 
(11. II) if and only if its Bloch wave expansion {Eg,Bg) G i^^(T^) satisfies G{Eg,Bg) = 
for almost all 9 G [0,lp with G from fl6.3p . A function u is in the invariant space of 
functions orthogonal in ^J^.o^V) these stationary solutions if and only if its expansion 
satisfies div(eo-E'e) = div(/io-B6») = m the sense of distributions. 



15 



6.3. Bloch spectral theory. Consider fixed periodic eo(a;), /io(x) G C°^(T^). The func- 
tion 

u{t,x) := e^* {E{x) , B{x)) 

satisfies MaxweU's equations (11. ip if and only if 

feoix) \fE{x)\ _ ( cni\\(E\ 
^^■^^ ^ V /io(a;) J V5(x) ) ~ \-cm\ V^y' " 

In the same way a 6'-periodic u = e^'^{E{x), B{x)) is a ^-periodic solution of Maxwell's 
equations if and only if _B is a solution of (16.51) and {E , B) is 6'-periodic. The function 
u is then a solution of (12. 2p . 

When each of eo and fio is a positive constant times the identity the change of variable 

' — ' ' — ' 1/2 1/2 

E, B = eQ E, j^Q B reduces to the case Cq = /Uq = I- For that case the eigenvalue problem 
is exactly solved in the next example. 

Example 6.6. In case eQ = fiQ = I the problem is translation invariant. Denote by the 
translation operator u{x) h-> u{x — i) acting on L^(Tg). The antis elf adjoint G commutes 
with Ti so the eigenspaces ofTg are invariant by G. 

For given k & 1? denote by E^ the subspace consisting of exponentials e*'^'^e*^'^(e, b) 
when e, b run in M? . Ek consists of eigenvectors of Ti with eigenvalue Choose 
the vector i so that the are rationally independent. Then the eigenvalues for distinct 

k are distinct, so the spectral decomposition ofT^ is L'^{Tq) = Q)±Ek. 

It follows that G{Ek) C E^ for all k & 1? . It suffices to diagonalize the restriction of 
G to each Ek. Compute 

r ( je.xik.x f^W _ ( curl\ / ie.^ / e'^-'^eX \ _ je.xjk.x ( i{k + e) M>\ 
\h) ) ~ V-curl yl V V^'^^'V ^ - e e q-^ ^ ^ j ■ 

To have eigenvalue X = iu it is necessary and sufficient that u is an eigenvalue of Gq G 
Hom(C'') 

One has eigenvalue if and only if both e and h are parallel to k + 6 This kernel has 
dimension equal to two. 

The orthogonal space has dimension 4 and consists of vectors with both e and h or- 
thogonal to k + e. Since {k + 9) J.h, {k + 9) A {k + 6) Ah = -\k + ^pb. Using (K^ 
compute 

-\k + e\'^h = {k + 6) A {k + 6) Ah = {k + 6) A (we) = -00% . 

Therefore for {k + 6) ^ there are two roots u = ±\k + 6\. Each has a two dimensional 
eigenspace generated by taking e ± {k + 9) and b = =p(A; + 6') A e. 

Return next to the general case. 
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Proposition 6.7. There is an infinite dimensional space of 9 -periodic solutions of (16. 5p 
with eigenvalue A = consisting of E and B with vanishing curls. The ^^{J^q)- 
orthogonal complement to this kernel satisfy div(eo-E') = div(/io-B) = 0. There is a constant 
C independent of 9 so that 9-periodic E, B satisfy 

\\E,B\\Hi(T'i^) < Cfllcurl E||i2(T3) + llcurl 5||^2(Tr3) 

(6.7) 

+ ||div(eo(a;)E)||i2(T3) + ||div (/io(x)5) 11^2(^3) + \\E, B\\L2(j3)j . 

For each 9 the X = iu for which (16. Sp has a nontrivial solution is a discrete set in 
M \ {0}. Each eigenspace is a finite dimensional space of smooth functions. The value 
is not an accumulation point of nonzero eigenvalues. The u are not bounded above and 
are not bounded below. 

Proof. The estimate is tlie key. Write u = e^^^u witli periodic u. Expand tlie latter in a 
Fourier series. The proof of the Lemma [3.31 yields (16.70 . 

Proposition 16.51 implies that the stationary solutions are curl free and their orthogonal 
complement is invariant under the flow by Maxwell's equations. Also that the complement 
consists of solutions satisfying div(eo(a;)-E') = div{fiQ{x)B) = 0. 

Decompose L^^^^(Tg) as a Maxwell invariant direct sum of stationary and dynamic 
states 

Ll,^oi^l) = KerG ©± T^dyn , Hdyn := {E,B : div(eo(a;)^) = div/io(x)5 = O}. 

The LI^^^^{Tq) unitary group e*"^ restricts to a unitary group on Hdyn whose anti selfadjoint 
generator is the restriction G\-n^^^. 

Estimate (16. 7p implies that (/ + G\-}i^^^)^^ is compact, hence has pure point spectrum 
tending to zero and total multiplicity in{|z|>(5>0} finite for all 5 > 0. Therefore the 
spectrum of Gj^^^^ is pure point and the total multiplicity in any bounded set is finite. 

Commuting with derivatives yields an inductive proof of an version of (16. 7p for 
1 < s G N, 

\\E,B\\hs+^tI) < G(s)( llcurl E||^s(T3) + llcurl fi||^.(Tr3) 
(6.8) ^ 

+ \\dw{eo{x)E)\\Hs^j3-j + ||div(/io(a;)5)||^s(T3) + \\E, B\\hs^j3) 

The smoothness of eigenfunctions for eigenvalues A 7^ follows. 

It remains to show that the spectrum is unbounded above and unbounded below. Define 
P to be the bounded strictly positive selfadjoint multiplication operator P{E, B) := 
{eE,fj.B). The eigenvalue equation is Gu = icuPu. It holds if and only if f = P^^'^u 
satisfies p-^/\G/i)P-^/^ V = ujv. Need to show that the eigenvalues are not bounded 
below and not bounded above. The u are bounded below (resp. above) by C G R if and 
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only if for all v so that P ^I'^v belongs to the domain of G, 

{P~^l''{GI%)P-^l''v,v) > C{v,v), (resp. <). 
This holds if and only if for all u belonging to the domain of G, 

{{G/t)u,u) > C{P^/\, P'^\) , (resp. <). 

This holds if and only if the spectral problem with = fiQ = I has eigenvalues bounded 
below (resp. above). Example (16. 6p shows by explicit computation that for = = I 
the spectrum is unbounded in both directions. Thus the ui are not bounded below and 
not bounded above. □ 

Proposition 16 . 71 implies that the discrete spectrum at fixed 9 consists of {0} with infinite 
multiplicity and a discrete set of possibly multiple eigenvalues \ = iu that we label 
according to their distance from the origin as in (11. 5p . Away from eigenvalue crossings, 
the functions 9 i— )■ Wj(^) are real analytic. Rellich's theorem shows that away from the 
crossings the associated spectral projections T\.j{9) are also real analytic in the sense that 
the unitary map of Proposition 16.41 intertwines them with an analytic family acting on 
Ll,m (see my 

For 9 fixed and an eigenvalue uj{9) ^ there is a finite dimensional space of eigen- 
functions e^^'^{E0{x), Bg{x)) G L^(Tg) and corresponding Bloch plane wave solutions 

of dEH) 

We assume the constant multiplicity Hypothesis II. 4[ 

From the analytic dependence of the operators it follows that the L^(Tg) orthogonal 
projection, n(^) G Hom(L^(Tg)) onto the nullspace of iu{9)AQ{y) — J2j ^j^j analytic 
and in particular of constant rank. 

7. The purely periodic case 

Fix 9 and a locally constant multiplicity eigenvalue u!{9) ^ 0. Denote by ¥Jyff) the 
kernel of L(ci;(^), ^, from Definitions 11.61 \i 9 ^ e^^'^il){x, 9) is a smooth function of 
^ on a neighborhood of 9_ with values in K(6') then il) is periodic with period 27r in x and 
smooth in its dependence on 9 for 9 ^ 9_. The function 

is a 6'-periodic solution of Maxwell's equation in the periodic medium e(a;),/i(x). 

Scaling the periodic structure to e{x / h) , fi{x / h) yields the corresponding rapidly oscil- 
latory Bloch plane waves 
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For a G C^(]R'^) and h « 1 the function a{{9 — 0_)/h) is supported in the domain of 
definition of u{9). Superposing nearby waves yields a Bloch wave packet for /i << 1 

i[o,ip \ h J 
Letting ( := {0 — 6)/h yields the exact solutions 

(7.1) = e'^-"/'* I ^f^, ^+/iC) e'*"(^+'^^)/'* e'"-^ a(C) c/C- 
7.1. The geometric optics time scale t ~ 1. 

Definition 7.1. Complementing Definition \1.5\ the corresponding transport operator is 
defined by 

I]){dt,d,) := dt + VA. 
The symbol is D(r, k) = t + V.k. 

The Taylor series in h of infinite and finite orders respectively are, 

(7.2) ^{yA + hO ~ i'{y.0) + Y.y 9,{yX). oo{e + hC) = uj{e) -VhC + h^k{hX) , 

where the sum on the right hand side of ~ is an asymptotic expansion as /i — )■ 0, not a 
convergent series. Then, 

('^•3) = e**-(^)/'^ e-^*^-« (l + Y^^htyHh. O) , 

where the last line uses a Taylor expansion of s t— )■ e**^*^'^'^-* about s = 0. Define 

v{x) := j e^^-^a(C) exc- 
use (D and O in dH]) to find 

(7.4) ~ e'^''^ w{h,t,x,x/h) , S := w(^)t + a;.^, 

and 

w{h,t,x,y) ~ Wo(i,a;,?/) + Wi(t,x,?/) H 

in the sense of Taylor series about h = 0. The leading term is 

woit,x,y) = v{x -Vt)i!{y,9) . 

The velocity is V = — V6»cj(^). It is not at all obvious from this definition that V does 
not exceed the speed of light. As our media are anisotropic, the speed may depend on 
direction. We recall the algorithm giving such anisotropic speeds (see [23], [19], [25]). 
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Denote by (t, ^) the dual variable of {t,y). The characteristic polynomial of Maxwell's 
equations is 

Its roots T{y,^) for ^ real define the characteristic variety. Define extreme roots Tmaxd/, 

by 

T-max(z/,0 := max{r : p{y,T,^) = O} . 

with a similar defintion for Tmin{y,0- The function Tmax is positive homogeneous of degree 
one in ^, and, 

(7.5) Tmax(Z/, -0 = -^min(l/,0- 

The set of velocities that do not exceed the speed of propagation is a convex set of vectors 
V whose support function equal T^a.x{y, —0; equivalently 

ngeM3\{o} {v : ^.v < r„,ax(y, -0}- 

In case of constant and isotropic permittivities rmax(l/, — = l^l/v^ yielding the 
classic formula for the speed of light. 
Define 

T"max(0 := max T^s,x{y,^). 
Then rmax(— is the largest speed limit for ^.v. 

Theorem 7.2. The group velocity V from Definition \1.5\ respects the maximum speed of 
propagation for each C, ^ 0, precisely for all E M'^, .^.V < Tmax(— 0- 

Proof. The Bloch spectral eigenvalue problem for periodic as opposed to 6'-periodic func- 
tions is 

{'t M.)) ^ ( ; ) - « 

where 11(6') is the projector (of constant rank) on the eigensubspace of the eigenvalue 
u{6). The constant multiplicity hypothesis implies that the eigenvalue and the spectral 
projector are analytic in a vicinity of 6_. Differentiating fl7.6p with respect to 6 in the 
direction of a covector ^ and multiplying on the left by 11(6') yields 

(7.7) ev«^»)nw('f) ^;^,)n(9) + n(9)(;^ 1'')n(e) = o. 

Take any eigenfunction E{y), B{y), normalized by 
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The quadratic form associated to (17. 7p yields 

The —T{y, ^) are the eigenvalues with respect to the postitive definite matrix diag {eo{y) , yUo(y)} 
so the min-max characterization implies that for each y 



-^A\(E{y)\ (E{y)\ ^ _^ (vniMy^ \ ( E{y)\ ( E{y)\ 
J {B{y) J ' {B{y) J ^ " ^ ^ ^o(z/) J \B{y) J ' [B{y) J ' 

- ^^(y) [B{y)h[B{y)P- 



Integrating and using f l7.8p and ( 17. 9p proves 

CVecoie) - r^in(0 > 0, 

which is the desired relation since V = —'V0Uj{6). □ 

An analogous result to Theorem 17.21 is proved in [3J , where a bound on the group 
velocity for scalar wave equations is given. 

7.2. The diffractive time scale t ~ 1/h. In (17. 3p . the expansion parameter is ht so 
when ht is not small, the approximation is not appropriate. For the diffractive scale 
/it ~ 1, one needs a refinement. Take the next term in the Taylor expansion in the 
exponent. Denote by q the symmetric quadratic expression 

(T.10) „c,c). tggcc,. 

t,j=l ■> 

Then, 

u{e_+eC) = uj{6) - hV.C + h'q{C,0/2 + /i^ /i^/3,(C) , 

i>o 

and. 

Introduce the slow time T = ht. The exact solution has the form 

Wih,r,x,y) := j ip{yA+hC) e^^^^^'^^/^ (.ihTj:,>,h^P,{Q g^-C . 
Taylor expansion in h yields 

(7.11) e'hTj:^>o'''^^^^^ = (^i + J2h^r^(rx))- 
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Using (17. 2p and (17. lip in the definition of W sliows tliat 

(7.12) W{h,T,x,y) ~ Yl h^^,{T,x,y), 

witli 

(7.13) wo{r,x,y) = ij{yA) J e'^"'^^'^^/' e''< a{0 dC . 
Tliis sliows that the solution has an asymptotic expansion of the form 

with W satisfying (TTT^ . 

Equation f l7.13p implies that wq is a tempered solutions (with values in K) of the 
Schrodinger equation 

(7.14) tdrwo - ^d'Mi){d.,d,)wo = 0. 

Though the function wq takes values in the finite dimensional space IK the equation f l7.14p 
is scalar. The constant rank hypothesis is crucial here. 

8. Block wave packets on a modulated background and t = C(l) 

This section considers solutions of the Maxwell's equations for times t = 0{1). This 
time scale is an essential first step in treating the diffractive case. In order that the asymp- 
totic description be nontrivial we allow lower order terms in the equations. In particular 
this includes the case of a possibly conducting medium with Ohm's law dissipation, 

(8.1) e{x,x/h) dtE = curli? — (r{t,x,x/h)E, ^{x,x/h) dtB = — curli?. 

Here a{t, x, y) is a nonnegative symmetric matrix valued function. The physics modelled 
is that where a ^ the medium reacts instantaneously to the field E by generating a 
current J = a E. Such an assumption is realistic only when the field E varies little on 
the time scales associated to the motion of electrons. The associated energy dissipation 
law is 

dt / {eE , E) + {fiB, B) dx = - {aE , E) dx < . 

The lower order tem is 0{1) and appreciably affects the fields for times t = 0{1). 
In this section only, we replace Hypothesis 11.11 by the following that allows (18. ip . The 
perturbations are larger by a factor than in the diffractive case. 

Hypothesis 8.1. T = C(l) hypothesis. 

3 

(8.2) = P^{t,x,dt,d^)u^ := dt{A^u^) + J^'^A,^^ + M'' . 
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The coefficients Aj, for j = 1,2,3, are the constant 6x6 matrices 

Jl 



(8.3) 
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The coefficient Aq and are of the form 

(8.5) A^it.x) = A'^o{x/h) + hAl{t,x,x/h) , = M{t,x,x/h) 

where Aq and Aq satisfy Hypothesis 



Remark 8.2. It follows that the growth rate C{m,h) from Theorem \4.1\ (and Remark 
is of order C(l) as 0. 



Motivated by the special case of purely periodic media in Section 17.41 and the linear 
case of Lax [22] (see also p5]) try the ansatz of two scale WKB type, 

(8.6) v\t,x) := e*^(*'^)/'^iy(/i,t,x,^), W{h,t,x,y) = WQ{t,x,y) + hw,{t,x,y) , 

where the Wj{t,x,y) are periodic functions of y with period 27r. The case when is a 
linear function of (t, x) is our principal interest 

s{t,x) = ujt + e.x, (w,^) gm^+^\{o}. 

For this phase the rays will be parallel straight lines and one finds Schrodinger type 
equations at the diffractive scale t = 0{l/h). 

Three identities are at the heart of checking the accuracy of the ansatz 

'iu 



(8.7) 

These yield 
(8.8) 

Then 
(8.9) 



curl 



/l|;e^^(*'^)/^iy(/i,t,a;,i/)] = e*^^*'^)/'^ + 9*) A'^qW {h,t,x,y) 



iS(t,x)/h w{h,t,x,y)\ = e*^(*'^')/'' ]-dyAW{h,t,x,y). 



—curl y 

h " 



P^{t,x,dt,d^)v'' = e'^^''^^/^ Z\t,x,x/h) 
(5. + f)A 



with 



A 







^ a, A 





h 



M 



W. 



Z'\t,x,y) = h^^r.i + ro + hvi + /iVs 



rj{t,x,y) . 
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Since one substitutes y = x/h, it would suffice to satisfy rj = on the subspace of (t, x, y) 
with X parallel to y. We achieve the more ambitious goal of choosing the Wj so that 
r_i = ro = everywhere. 

8.1. The leading order term. The leading two orders in are 

+ ro = h'^ l.{u,e,y,dy)W + {M{uj,e,y,dt,d^,dy) + iuAl+ M^W , 
where L is from (11. 6p and 

(8.10) M{y,dt,d,) := A',{y)dt - 

The highest order term in (18. 9 p is 

(8.11) r_i = h{u,e,y,dy)wo. 

In order that r_i = have nontrivial solutions, it is necessary that 

keihiu, e,y,dy) ^ {0}. 

According to the Floquet-Bloch theory of Section El L(a;, 6, y, dy) has a nontrivial kernel 
on periodic functions if and only iuj = ioj{6) is an eigenvalue of (II. 4p . From now on we 
make the choice of d_ and u}{6) so that the constant multiplicity hypothesis 11.41 is satisfied. 

Definition 8.3. In addition to DefinitionlTE denote by Q e Hom(if"(T3); H'~^\Tl)) the 
partial inverse of L defined by 

Qn = nQ = o, QL = LQ = /-n. 

The equation r_i = is equivalent to G K := kerL, that is 

(8.12) Uwo = Wo. 
Using the definitions of L and M, the term ro is given by, 

(8.13) ro = hwi + (M+ iuAl + M)wo, 
so tq = if and only if, 

(8.14) hwi +(M+ iuAl + M)wo = 0. 

Equation (I8.14p involves both wq and Wi. Since L is selfadjoint on L^(T^) its range is 
perpendicular to its kernel so IIL = 0. This is true only because 11 is the orthogonal 
projection on the kernel IK for the L^(T^) scalar product (not for the other scalar product 
in Definition II. 6p . The equation Ilro = yields an equation for wq alone, 

n(M + iujAl+ M)wq = 0. 

Taking into account (I8.12p . this is equivalent to 

(8.15) n(M + iuAl+ M)Uwo = 0. 



f d^A\ 
V -d.A J ■ 
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Proposition 8.4. For any w(t,x,y) G C°° , 

(8.16) nMHw; = U A°U {dt + VA) w , 

with the group velocity V from Definition } 1.5[ The operator TlA^Tl is a linear isomorphism 
ofK. to itself. 

Proof. Prove the last sentence first. Since K is finite dimensional, it suffices to prove 
injectivity. If 11 11 A; = then 

= {UAoUk,k) = {AoUk,Uk) > c\\Ukf , c>0, 

since Aq is strictly postive. Therefore k = Uk = proving injectivity. 
From the definition of 11 and M one automatically has for arbitrary w, 

d 



with matrices aj{y). It suffices to compute the aj. This is done by computing the dif- 
ferential operator on the test functions tip, and Xjip, with G K. Applying to tip and 
setting t = yields 

(8.17) aoU = nA°n . 

Applying to Xjip and setting Xj = yields 

(8.18) n = — n ^ ^ n , {ei, 62, 63} is the standard basis of . 

The identification of aj from fl8.18p and (18.171) requires first order perturbation theory as 
in (18.191) of the next proposition. Second order perturbation theory as in (I8.20p is needed 
for diffractive geometric optics. The identites are proved by differentiating the identities 
IIL = and LII = 0. We refer the reader to [2], [25] for detailed proof. 



Proposition 8.5. Suppose that 9 and cu satisfy the constant multiplicity hypothesis iLj 



Suppose that the coefficient A^ and 9 depend smoothly on a parameter a with their unper- 
turbed values attained at a. With ' denoting differentiation with respect to a, the following 
perturbation formulas hold, 

(8.19) HL'n = 0, 
and 

(8.20) nL"n - 2nL'gL'n = o. 

Returning to the formula for a^, use (I8.19|) with a equal to the j^^ component of 9 and 
' = d/d9j. Then, 



25 

The above identity in combination with f l8.19p and flS.lSp yeilds 

a,n = -nA°(y)n|^. 

This together with fl8.17p completes the proof of Proposition 18.41 □ 
Recall Definition 11.71 The map 7 inherits the regularity of the coefficients, 

(8.21) a°^7 G L°° ; HomK) . 

Then r_i = Ilro = exactly when wq = Hwq satisfies the transport equation 

(8.22) (^dt + VA + l{t,x)'^wo = 0. 
The equation ro = is equivalent to the pair 

nro = 0, Qro = 0. 
Equation (18.130 shows that Q tq = if and only if 

(8.23) {I-U)wi = -Q(M+ iujAl + M)wo. 
The choice of Uwi does not infiuence ri,rQ. Choose 

(8.24) nwi = 0. 

Theorem 8.6. If g e C'=^{R^;K) there is a unique Wq G C°°(M; //""(M^; K)) satisfying 
flS:^ with wo{0) = g. Define Wi and by ([H23D, flOD and ([HS]) respectively. If is 
the exact solution of P'^u'^ = with m^L ^ = v'^l, n, ^^en for all a 

•' lt=0 lt=0' •' 

sup \\{hdt,.r{u^-v'')\\^,.^,. < C{a)h, 0<h<l. 

te[o,T] ^ ' 

Proof. I. Estimate for P^v'^. Use (18. 8p together with (18.91) and the fact that the equa- 
tions satisified by the functions Wj guarantee that r_i = tq = 0. Therefore for h g]0, 1[ 

ll'^M,j/^''ILoo([o,r]xK3xT3) ^ C{a)h, supp Z'' C {{t, x + Vt,y) : x e supp ^} . 
This implies the fundamental residual estimate 

||(/^9,.)"(pV)||^^(j^_^,^^3) < C{a)h, he]o,i[. 

Using the compact support one has 

(8.25) \\ihdt,.r{P'v')h^^^o,Ty,LHm ^ ^'(a) /i, /iG]0,1[. 

II. Stability for P^. For T > and m fixed Theorem 14.11 shows that there is a 
constant C = C{T, m) so that for < t < T and h g]0, 1] 

(8-26) |h(t)|L,, < C (|h(0)|L_, + £ \\P'w{s)\\^^^ ds) . 
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III. Combining. Apply f l8.26p to := u'^ — to find 

ll(-'-^'')WlL,. < ^(ll(-'--')(o)ll™,. + I - -'MU ds 

= C(\\{u'-v'm\\, + f\\P'v\s)\\ds 



^(ll(-'--')(o)IL„,. + 0{h)) 



where the last step uses the residual estimate fl8.25p . 

It remains to show that = — satisfies for for all m, 

(8.27) lk'(0)|U,^ = 0{h), 

To do that it suffices to show that for all < A; G N and s, 

(8.28) \\ihdt)''w\0)\\H^^^u^) = 0{h), h^O. 

Since the initial values vanish, the case A; = is automatic. Prove f l8.28p by induction on 
k. Suppose (I8.28P known for < k < k and all s. We prove it for k + 1 and all s. Begin 
with the identity 

j 

The first term on the right is hdtw^. Therefore 



The inductive hypothesis implies that \\{hdt)-'^^u(fi)\\Hi{9:i) = 0{h). This completes the 
inductive proof of fl8.28p and therefore the proof of the theorem. □ 

9. Block wave packets on a modulated background and t = 0{l/h) 

This section considers solutions of the Maxwell's equations with coefficients and lower 
order terms satisfying Hypothesis 11.11 so as to become pertinent at t ~ l/Zi. 

Remark 9.1. With permittivities satisfying Hypothesis \l.l\ the growth rate C{m, h) from 
Theorem |^.i| satisfies C{m,h) < c{m)h with a constant c{m) independent of h g]0,/;,[. 
The time evolution is uniformly hounded so long as the product t x h remains bounded. 



Motivated by the special case of purely periodic media in Section 17. 4[ the ansatz ex- 
pected to be valid for times t = 0{l/h) is of two scale WKB type, 

(9.1) v'^it.x) := e'^'^''''^/''w(h,ht,t,x,^^ , S{t,x) = ut + e.x, 

W{h,T,t,x,y) = wo{T,t,x,y) + hwi{T,t,x,y) + h'^W2{T,t,x,y), 
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where the Wj{T,t,x,y) are periodic functions of y with period 2it. To go further in time 
requires the additional corrector W2- In order to preserve the relative ordering of the size 
of the terms for t = 0{l/h) it is crucial to insist that the Wj{T, t, x, y) are suhlinear in t, 

Wj{T,t,x,y) 



lim 







uniformly in T, x, y. We construct profiles satisfying a stronger hypothesis. 
Compute using fl8.7p 



(9.2) 



(9.3) 



Then 



P''{t,x,dt,d,)v^ = e'^^''''^/^Z\hr,t,x,x/h) 



with 



Z\r,t,x,y) :-- 



{d^ + i9/h)A 

{d^ + i9/h)A 







idyA 
-\dyA 



hM 



W . 



(9.4) Z'^{t,x,y) = h ^ r_i + Tq + hri + h"^ r2 + , rj = rj{t,x,y). 

As in the case of t = 0{1), we achieve the ambitious goal of choosing the wj so that the 
leading r_i = tq = ri = everywhere. This reduces the residual to 0{h'^) which allows 
us to justify for times t = 0{l/h). 
One has 

r_i = Lwq and tq = hwi + Mwq 
with L and M defined by (11.61) and (18.101) respectively. The next coefficient is 

{ie + dy)A' 











^-{^e + ^y)A ■ • J^^ + Mwo 
= hw2 + Mwi + Nwo , 

where 

(9.5) N := 9rA° + zwAj + M. 

As in Section ^ r_i = and 11 Tq = if and only if 

riwo = , and {dt + V.dr,)wo = , 
with V as in Definition II. 5[ Thus there is a reduced K valued profile Wq such that 

wo{T,t,x,y) = wo{T,x -Vt,y) . 
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In order to determine Wq one needs a dynamic equation in T. The equation Q tq = 
yields 

(9.6) (/ -U)wi = -QMwo, 
and thus 

(9.7) {dt + VA){l -Il)wi = 0. 

The equation Ilri = yields the Schrodinger equation determining the dynamics of wo- 
Plugging (19. 6p into the equation Ilri = yields 

(9.8) n M n^i - n M Q M n^o + n n Uwq = o . 

Proposition 18.41 and (19.71) imply that equation (19. 8p is equivalent to 

(9.9) (HNn-nMQM U)wo = -U A° U{dt + VA)Tlwi . 
The next proposition identifies the operator on the left hand side of (19. 9p . 

Proposition 9.2. On smooth functions w{T,t,x,y) that satisfy {dt + V.dx)w = 0, 

(9.10) {UNU - nMQMn)u; = UA^Quldr + d^uj{d^,d^) + -f{t,x)]w . 



2 

Remark 9.3. It is surprising to find that the operator 

(nA°n)-i(nNn - nMQMn^ 

acting on K valued functions has leading terms that are scalar. Coupling only occurs 
through the zero order term 7. A related zero order coupling occurs in §6 of [4j 

Proof. With k = (fci, k2, k^) G fixed, 6 := 6 + ak, a G M, differentiate L(a;, 6, y, k) with 
respect to a to find 

= ^kAl^eu - '^^^ = zM{y, d„ k) - lA^ ©(9,, k) , 

L" = iikdefcoA^Q . 
Using the first identity yields 

n L' Q L' n = -n M Q M n - n A° ©(Si, k) q a° ©(9t, k) n . 

Applying (18120]) yields 

(9.12) ^ n L" n = -n M Q M n - n a° ©(St, k) q a° B{dt, k) n . 

Equation (19.120 and the definition of N in (19.50 , give 

'hnh - nMQMn)w7 



(9.13) 



n 



drAl + iuoAl + M+-i {kdefuoAl + A° ©(9i, k) Q A° 0(9*, k) 



Yiw. 
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Next replace k by and simplify the right hand side of 09.131) using 
This yields 



□ 



9.1. Ray averages. In general it is impossible to satisfy f l9.9p exactly since all the terms 
are annihilated by dt + Vdx except the 7(t, x) term. If the coefficients are constant on 
group lines the 7 term is too and one can construct infinitely accurate solutions (see [2]). 
In the present case we replace 7 by its average on rays to find a solvable equation. Then 
estimate the error. For that estimate we impose an assumption on 7 slightly stronger 
than the existence of ray averages. The material is recalled from [2] where the proofs and 
additional information can be found. 

Assume that 7 e C~(Mi+3 ; Hom(]K)) satisfies f lST^T]) and that the averages on rays exist 
as in (11. 8p . It follows that the 7 is smooth and that the ray averages of the derivatives of 
7 exist uniformly on compacts and satisfy 



(9.14) 



lim 



T 



- / &>d^,^{t,x + Vt)dt 



X 



0. 



We need more than this. The function 7(x) is the average on the ray intersecting t = 
at X. The ray passing through the point {t, x) intersects t = at x — Vt. The function 
which assigns to (t, x) the average value of 7 on the ray through (t, x) is equal to ^{x — Vt). 
The function that subtracts from 7(^,0;) its average on the group line through (t,x) is 
equal to 7(t, x) — 7(0; — Vt). 

Consider the Hom(K) valued solution g of the transport equation 



(9.15) 
Then 



+ V.d^^g = 7(t,x)-7(a;-Vt), 9^^0 = 0. 

g{t,x) 



'y{s,x + Vs) ds — 7(x), 



X := X — Vt . 



Assumption (II. 8p is equivalent to the fact that this is o(l) as t — )■ +00, 
(9.16) lim sup = 0. 



Assume that 7 satisfies the ray average hypothesis in Definition 11.8 
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Example 9.4. i. If 'y{t,x) = f{i{t,x)) where f{6) is a smooth periodic function of 
arbitrary period and i is a linear functional then the ray average hypothesis is satisfied 
with (3 = 0. 

ii. If Ai : R-'^^^ — 7- is linear and satisfies the (generic) small divisor hypothesis 

3C>0, mGN, VnGN^^ n.M{l,V) ^ \{n.M{l,V)\ > C {nl^"" , 

then, for h{9i, ... ,9m) € C°°(T^^) the quasiperiodic function 'y(t,x) = h{M.{t,x)) satisfies 
the hypothesis with /3 = (see [17j ). 

iii. Consider smooth almost periodic 7 of the form 

(9.17) ^{t,x) = 

where a^i vanish for all but a countable family of rj and satisfy 

(9.18) VnGN, 5^ (^)" Isl < 00 [t]) ■= {1 + \r]\'^y/\ 

V 

Then •yit, x) — 7(x — Vt) = y^_^Q ^^■n.{t,x) _ Suppose that there is an a > so that 
for all n 

(9.19) (^)" h'^l = ^(^")' '^^0- 

0<|r?.(l,V)|<<5 

Then the ray average hypothesis of Definition \1.8{ holds with (3 = a/{a + 1). 
9.2. Using ray averages. Rewrite (19.91) as 

(9.20) (dr+^t dlu{d,,d,) + j{x-Vt)^wo = -{dt + V.d,)Uw,- {y{x)-^{x-Vt))wo . 

This equation is satisfied by choosing wq and Uwi so that both sides are identically zero. 
The left yields equation (11.91) . The initial value, u;o(0,x) G S{M.^ ; K) is arbitrary. 

Lemma 9.5. For any f G 5(M'^; K) there is a unique solution C(7~, x) G C°°(M7- ; iS(]R^; K)) 
to the Schrddinger initial value problem 

(9.21) (^OT + lidMd.,d,)+j{x)y = 0, C(0) = /■ 

For each m,s,r G N, \a\ < m, \k,\ < s there exist constants c{m, s,r),b{m, s) so that for 

aiir 

\\x^d:dJrar)\\LHR^) < c(m,s,r)e(^+'('"'^)"^"M'-+-)^ (1 + imwr^o^Y ^ . 

KI=o 

Proof. We present only the a priori estimates. Multiplying the equation by ( and taking 
the real part yields 

[ \C\'dx + ^ [ {^CX)dx = 0. 



Estimate 

^ I {lC,Odx< II7IL00 / Kl'dx, 



lan' dx<^ I \f\Ux + 117IL. / icmr dx. 

Jm.^ 



and integrate with respect to T to find 

If \ar)\'dx<l 

Gronwall's inequality gives 

(9.22) llC(r)|U2(^3) < ell^ll.-^ II/IU. 

The Duhamel relation 

d^^Cm = ir(r,0)9,,C(0) + r K{T,s)P{d^^C){s)ds 

Jo 

= KiT,0)d,J + f KiT,s)i-d,^^)as)ds 
Jo 

then yields the estimate 

||5x,C(r)|U2(i,3) < e'l^".-^ ||9.^./|U2(M3) + rell^ll^-(^-^)||V7|U°ci|C(s)|U^(M3)rfs 

^0 

This combined with (19.221) gives 

||C(r)bi(M3) < 2 ell^ll^°°^ ||/||^,i(M3) + rell^"^-(^-^)||V7|U-||C(5)l|L^(M3)rfs. 

Jo 

Now apply Gronwall's inequality to get 

IIC(r)ii^,i(M3)<2e"^iw -^11/11^,(^,3). 

By induction one proves that 

sup II (9.)" cm 11^,(^3) < 2e''(™)IRII— ^l|/l|^,.(M3). 

\a\<m 

Let us prove the weighted estimate x'^d^. The commutator of Xj with the Schrodin 
operator S in (19.211) is the first order scalar differential operator 

[S,Xj] = i^{dluj)ijda,^ . 
I 

Therefore, for |a| < m 

xjd^ar) = K{r,o)x,d:m + f K(r,s)s{x,dP^c){s)ds 

Jo 

= K{T,0)x,d:f + [ K{T,s)[S,x,]{d:Ois)ds, 

Jo 

which yields 
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By induction one proves, for \k\ < s, 

(9.23) \\x''d:ar)\\L^m < c(m,s)e(^+'("^'^)"^"M.™+».oo)^ ^ 11x^11^'-+^- • 

K|=o 

The time derivative commutes with the Schrodinger operator, therefore, for each r G N, 
satisfies the same equation as ( with initial condition 

\dKm<cir)ii + \mwr,^rJ2\^rf\- 

j<2r 

Apply flOSi) to 9f C to find 

||x«a:9fC(r)|U^(M3) < c(m, s,r)eM"^'^)"^"w-+»-)^ (1 + ||7|k^.o.)^ ^ \\x' fy^^s^.^.r . 

\e\=o 

□ 

Given Wo{0) = / G 5(M'^; K) choose Wq the solution provided by Lemma 1^31 It satisfies 
for all a, 

(9.24) {x,dr,t,.rwo G L°°([0,T]r x x ; K) . 

Setting the right hand side of fl9.20l) equal to zero yields an equation that is solved using 
a Hom(K) valued integrating factor g{t, x), 

(9.25) Ilwi= git,x)wo, 

where g is the solution of (19.151) . The ray average hypothesis with parameter < /3 < 1 
yields estimates for the derivatives of g and therefore those of Hwi, 

{ty^ (x , dt,.riUwi) G L°°([0, T] X Mi X R^. ; K) . 

The component (/ — Il)wi is given by (19.61) in terms of Wq so (19.241) implies, 

{x,dr,t,.,yr{I-Il)wi G L°°([0,r]xMiXM3 xR3), 

with wi is periodic in y. This completes the determination of wq and wi. At this stage 
one has r_i = = Ilri = 0. We choose W2 to that (/ — n)ri = 0. The latter equation 
holds if and only if 

(9.26) {I-U)w2 = QMwi + QNwo- 

This determines (/ — 11)^2. On the other hand, Ilw2 does not affect the profiles r_i, tq, ri. 
It is chosen equal to zero, 

(9.27) Uw2 = 0. 
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The estimates for Wo,Wi imply that the y-periodic W2 satisfies estimates analogous to 
those of wi so, 

(9.28) {t)-^ {x , dr,t,.,yrw, G L°°([0, T] x x x R^), j = 1, 2 . 
This completes the determination of the profiles so that 

(9.29) r_i = ro = ri = 0. 

Theorem 9.6. If f e S{R^;K) there is a unique wq G C°°(R; 5(M^ K)) satisfying ^Ml 
with Wo{0) = f. Define Wi, W2 and by ([nSD, dHSSD, and dH]) respec- 

tively. If is the exact solution of P^u^ = with u'^\^_q = v'^\^_Q, then for all a 

(9.30) sup \\{x,hdt,.r{u^-v'')\\ < C{a)h^-^ < /i< 1 . 
te[o,T//i] ^ ' 

Proof. Let m G N. The bound fl9.28p and the identity (19.291) yield the residual estimate 

(9.31) ||(t)-^(a;,/.9,.)"(PV)||^^(j„^^/,j^^,) < C{a)h\ he%l[. 

This combined with Theorem 14.11 and Remark 14.31 shows that there exists a constant 
C(m,T) such that for all t G [0,T//i] 

II(«'--')WIL,. < ^K^)(ll(«'--')(o)IL,. + I ll^'-'(^)IL,.^^) 

(9.32) = C(m,T,/3)(||(«^-t;'^)(0)|L^ + 0{h'-P)) . 

For the first term in the rhs of (19.321) we follow the proof of Theorem 18.61 which yields in 
this case 

ll(^'--')(o)IU<^(^'"')- 

The error estimate with polynomial weights requires an additional weighted stability 
estimate. We will use the following notation. For a function u{t,x), integers i,m > 0, 
define 

\Ht)\\i^^^h ■■= h'^ih9t,xTu{t)\\^^^^^^. 

\l3\<e,\a\<m 

Lemma 9.7. Under the assumptions of Theorem \4.1\ there exist constants c{i, m), C{m, h) 
such that 

(9.33) \Ht)h,m,h < c(£,m)e^(^'™)^('"''^)* |k(0)||£,m,/^. 

Proof. It is a commutator argument resembling the proof of Theorem 14.11 The proof is 
inductive in i. The case ^ = is provided by Theorem 14.11 Assume that fl9.33p holds for 
£ = i. The case £ + 1 is proved by applying the inductive hypothesis to XjU for j = 1, 2, 3 
and using the commutator relation 

[P\x,] = A, J = 1,2,3. 
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An application of Gronwall's inequality finishes the proof of the lemma. □ 
Using the weighted estimates as above proves the theorem. □ 
Proof of Theorem 11.91 Write 

The preceding theorem estimates the first term. It suffices to show that the differences 
— and — have similar upper bounds. The first follows from the stability for 
Maxwell's equations proved in Lemma 19.71 The second follows from Lemma 19.51 □ 
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